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Abstract 

After reviewing the main results obtained within a model for the in- 
tersection of two perpendicular flows of pedestrians, we present a new 
finding: the changeover of the jamming transition from continuous to 
first order when the size of the intersection area increases. 
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1 Introduction 




Figure 1: Intersection of two one-way streets of width M. The blue particles 
(►) move eastward and the orange particles ( ) northward. The parameter a 
determines the particle injection rate. The region bordered by the heavy solid line 
is the 'intersection square'. Figure taken from |13] . 

In this talk we will deal with crossing flows of pedestrians, modeled as 
hard core particles that move on a lattice. The geometry of interest to us 
is shown in Fig.[Tl it represents two intersecting streets of width M and, in 
principle, infinite length. There are two kinds of particles, those moving east 
(blue) and those moving north (orange or red); each move covers a single 
lattice distance. When an eastbound and a northbound particle have the 
same target site, there is a question of which one has priority. The answer 
is given by the rules of motion of the model, that is, by the particle update 
algorithm. 

Frequently used algorithms are random sequential update, parallel up- 
date, alternating parallel update, sublattice update, random shuffle update, 
and so on (some definitions are given in [1] and in [2]). Recently we were 
looking for an algorithm that would (i) make every particle advance, as long 
as it is not blocked, at unit speed, and (ii) provide a natural answer to the pri- 
ority question. Wishing to avoid certain disadvantages of existing algorithms 
we were led to introduce the frozen shuffle update [H |3]. While operating in 
continuous time 

• this algorithm assigns to each particle i injected into the system a phase 
Ti G [0, 1] which is the fractional part of the time ti of injection; the time 
intervals during which an injection site is empty {i.e., between the departure 
of an occupying particle and the arrival of the next one) are i.i.d. exponential 
variables of average 1/a ; the injection probability a = 1 — e""^ during a unit 
time interval is a convenient control parameter; 
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Figure 2: (a) A single lane of finite length L with an injection probability a 
and an exit probability (3. (b) Case M = 1 of Fig.[T} two intersecting lanes. 
The exit probability from the intersection site is unity; however, their mutual 
obstruction leads for each lane to an effective transit probability /3| = ^ through 
the intersection site. 

• during the nth unit time interval, n < t < n + 1, this algorithm visits 
the particles in order of increasing phases and advances particle i at time 
t = n + Ti at the condition that its target site be empty. 

This talk is centered around the crossing street system with frozen shuffle 
update; we will nevertheless stress, where applicable, the robustness of our 
results under change of algorithm. 

2 One-dimensional model 

On the one- dimensional lattice of Fig. [2^ the model described above reduces 
to a totally asymmetric simple exclusion process (TASEP). We label the sites 
X = —L, —L + 1, . . . , —2, —1, the origin being chosen at the exit. For an exit 
probability (3 = 1, due to the specific way [3] of injecting the particles, none 
of them ever blocks its successor; all particles will then traverse this lattice 
at unit speed, ffrcc = 1, and will be said to be in a state of 'free flow'. The 
statistics of the free flow is fully known. In particular, the particle density 
Pfree and the particle current Jfree = "ffreePfree are given by 

Pfree(a) = ^frcc(«) = . (2.1) 

1 + a 

For an exit probability /3 < 1, whenever a particle wishes to move off the last 
site and thereby exit the system, this step is executed only with probability 
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Figure 3: Intersection of two single lanes, (a) A waiting line of length at 
the intersection site divides the horizontal lane into a free flow and a jammed flow 
domain, (b) The random phases of the particles in the horizontal lane define a 
division of them into platoons. A heavy vertical red bar has been placed to the 
left of the last particle of each platoon. Those in the jammed flow domain are 
compact. 
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Figure 4: Monte Carlo and analytic result (|3.2p for the reflection coefficient R = R 
coefficient for M = 1. Figure adapted from [9]. 
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/3. If refused, the particle stays where it is and waits a unit time interval until 
it can make its next attempt. The blocking of the exit may lead to blockings 
further down the lane and may create an intermittent or permanent waiting 
line. 

A particle undergoing a blocking by its predecessor will be said to belong 
from that moment on to the waiting line. Let us denote the site of the 
leftmost particle in the waiting line by a; = —i{t), where i{t) = 1, 2, 3, . . . , L; 
we set i{t) = if none of the particles in the system has ever been blocked. 
The site x = —i{t) separates a 'free flow domain' to its left from a 'jammed 
domain' to its right; we will therefore say that it is the location of a domain 
wall. In the moving frame that has the domain wall as its origin, the density 
profile to the left is strictly constant and equal to Pfree, whereas to the right 
it decays after some weak oscillations rapidly to a higher value Pjam- A 
study of this profile was carried out in Ref. [3] and has, incidentally, raised 
an interesting question of the validity of domain wall theory [3 El E] for 
this problem. 

The domain wall i{t) performs a random walk. At fixed /3 we expect that 
for low enough a it is localized within some finite penetration depth $,{a) 
from the exit (the system is in a state of free flow), and that for high enough 
a it is localized within ^(a) from the entrance (the system is jammed). In 
the vicinity of a critical value a = adP) the penetration depth becomes of 
order L and in the limit L — )■ cx) a sharp critical point Q!c(/3) arises on the a 
axis. 

Whereas the particle density pfree and the current Jfree in the free flow 
domain are known, the analogous quantities pjam and Jjam in the jammed 
domain have to be calculated. The key to the exact solution resides in the 
concept of a platoon, defined as a maximal sequence of successive particles 
having increasing phases (see Fig.lSb). This concept is therefore linked to the 
frozen shuffle update algorithm. A platoon is said to be compact^ if, at any 
integer instant of time t = n, there are no empty sites between its constituent 
particles. An elementary calculation [3] shows that the average number u of 
particles in a platoon is given by 

i = l + l-i. (2.2) 

u a a 

The argument leading to pjam and Jjam begins by considering the exit at 
X = in contact with the jammed domain. It may be shown [3] that in the 
jammed domain 

• each platoon is compact; 

• two platoons are separated by either or 1 empty site and the density 
of the empty sites is /3/(i^ + /3) = 1 — pjam- 

^In our earlier work we reserved the name 'platoon' for what we now call 'compact 
platoons'. 
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Both features are illustrated in Fig.[3]3. Employing these properties we 
deduce that Jjam = /Spjam with the current Jjam given by the equation 

1111 , , 

(2.3) 



Jjam Jfrcc P ^ 

The current J effectively passing through the system is J = min( Jfrec, </jam) 
and hence Eq. (12.31) shows that the critical jamming point occurs for 

ac = /3. (2.4) 

Although Eqs. (12. 3 p and (12. 4p are elegant and simple, we have found no quick 
way to see that they must be true: the system does not have the particle-hole 
symmetry that facilitates the analysis of certain other TASEPs. 

A final remark is that, unlike with other updates, the current Jjam here 
depends not only on /9 but also on a, which comes in through the platoon 
structure. For the intersecting streets to be studied below we have /9 = 1, 
but Jjam will continue to depend on a. 



3 Two crossing lanes: the case M = 1 

We now consider the crossing of the two single lanes shown in Fig. 12b, where 
the exit probability from the intersection site is unity; this is Fig. [T] for the 
special case M = 1. By an 'exact solution' we will mean an exact expression 
[such as (12. 3p and (12. 4p ] for the critical point and current J as a function 
of the injection probability a. We expect again that there is a critical value 
a = such that J[a) = Jfrcc(«) for a < ttc and 

J(«) = Jjam(«) 

= [1 - R(a)]Jfroc(a), a>ac, (3.1) 

and the challenge is to calculate and Jjam(ct) in the jammed phase. The 
second line of (13. ip is a rewriting which shows that Jfree and RJfrec are anal- 
ogous to an incident and a reflected wave, respectively and defines the re- 
flection coefficient R; in Fig.|3^ the reflected wave has advanced to the point 
marked i{t). Since R(a) vanishes for a < ac axis and varies between and 
1 for a > ac, we may employ it as the order parameter of the transition. 

Fig.H] shows the behavior of R(a) for the system of two crossing lanes of 
Fig.[2J This case is in fact exactly soluble [S]. It leads to c^c = | and 



^ / ^ 2tt ± 1 , , 

f^^«) = ^5;7-T^;— ' «>«c = i (3.2) 



Comparison of this value o^c = | to Eq. (12.40 shows that it is as if each 
lane exerts a blocking effect on the other one equivalent to an effective exit 
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probability that we will call /3j (for a reason to become clear) and that takes 
the value /3l = ^. The curve fl3.2p is shown in Fig. HI 

It is worthwhile to note that in the jammed phase the front of the reflected 
"wave" propagates at a constant average speed, 

{i{t)) = VRt, a>a,, (3.3) 

in which (. . .) is an ensemble average and where VR{a) and R(a) are related 

by m 

= aR + (!-«). • ^'-'^ 

To conclude we remark that the same system of two crossing lanes can 
still be solved exactly |8] for unequal entrance probabilities ai and a2 and 
exit probabilities Pi and P2 (real ones, not the effective one mentioned above) 
less than unity. 



4 Streets of width M > 1 
4.1 Theory 

For crossing streets of width M > 1 we have not so far found any exact 
solutions. The division of the sequence of particle into platoons seems no 
longer helpful. We do not exclude that somebody can find the solution for 
a 2 X 2 intersection square, or perhaps for the asymmetric case of an M x 1 
square. However, if the general case seems out of reach theoretically, accurate 
numerical studies are possible that reveal interesting properties. We prepare 
the ground by defining the order parameters that will be relevant. 

We will restrict ourselves to two crossing streets of the same width M. 
We will number the lanes from the inner ones outward by an index m = 
1,2,..., M. To each m there corresponds an eastward and a northward lane 
and by symmetry the two are statistically identical (we never found any hint 
of symmetry breaking). In the mth lane we now have a relation analogous to 
(13.11) but augmented with the indices M and m indicating the street width 
and the lane number, 

J^(a) = [l -<(«)] Jfreo(«), (4.1) 

and a similar generalization of Eq. (13.41) . The are now the order param- 
eters. We may formally write the reflection coefficients as [9] 



which for M = m = 1 and Pi = ^ reduces to (13. 2p . and where the P^ 
now has the interpretation of an effective exit probability from the mth lane. 
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However, for M > 2 we can find the (or equivalently the (3^) only by 
simulation. 



4.2 Simulation: memory boundary conditions 

Simulations are necessarily carried out on finite lattices. A finite value of 
L (see Fig.[T]) sets an upper limit to the length i{t) of the waiting line and 
will cause a rounding of the transition. We have conceived |9] an algorithm 
allowing the simulation of systems having L = oo, using only a finite number 
of variables. The trick is to consider the M x M interaction square with 
special boundary conditions, termed 'memory boundary conditions'. For 
each lane we keep track of the particle positions in the intersection square 
plus one extra variable which, essentially, is the length of the waiting line in 
that lane. This eliminates all finite length effects and the sharpness of the 
transition point increases with the duration of the simulation. 



4.3 Jamming for street widths M < 24 



We carried out simulations of the intersecting streets of widths up to M = 24. 
It appears [H] , at least for the values M < 24 that were investigated, that the 
mth lane undergoes jamming at a critical point and that ajj < a^f_i < 
. . . < af^. We call the principal critical point. Fig.|5l obtained with 
the traditional finite L simulation algorithm, shows a snapshot of the system 
with M = 10 and L = 15 which has its inner lanes jammed and its outer 
ones in a state of free flow. Fig. |6] shows the ten reflection coefficients [a) 
for the M = 10 system, which by (14. ip are directly equivalent to the currents 
J^{a) [this one and all further simulations were carried out with the memory 
boundary conditions of section i.e., for L = oo]. 



It appears that the principal critical point a^!^ decreases with M and that 



its behavior is very well approximated by 

as shown in Fig.[71 where A = 1.287 and B = 2.306. The high precision of 
these results is due to the elimination of finite length effects. 



4.4 Jamming for larger street widths 

The question of the M — oo limit of the critical point was asked also in the 
context of the BML model [TOl El] , but none of the authors has been able to 
state whether or not this point goes to zero in that limit. On the basis of our 
above results one might guess that (14. 3 p is the correct asymptotic law and 
hence that for the model of this work the principal critical point does tend 
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Figure 5: Snapshot of intersecting streets of width M = 10. The inner lanes are 
jammed and the outer ones in a state of free flow. Figure taken from [9]. 




Figure 6: Reflection coefficients for M = 10. Figure taken from [9]. 




Figure 7: Principal critical point as a function of M. Figure taken from [9]. 
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to zero with increasing M. However, when the simulations are pushed to 
larger lattice sizes, a novelty appears [12]. This new phenomenon is already 
suggested by the fact that the initial slope of R^{a) is of increasing steepness 
as M becomes larger, and that infinite steepness would correspond to a first 
order transition. Hence for growing M the transition seems on its way of 
becoming first-order. This is confirmed by further investigation. In fact, 
for large enough M and when a increases, the free flow phase appears to 
become metastable: by means of a nucleation mechanism it may irreversibly 
turn into a jammed phase. This is exemplified in Fig.[8l A system of linear 
size M = 100 is started at time t = in a free flow configuration with 
a = 0.08729, both free flows having just arrived at the entrance of the empty 
interaction square. After a transient of no more than a few hundred time 
steps the system settles in what seems to be a stationary state in which it 
stays for the first 800 000 time steps. Fig.[8K shows the particle configuration 
on the intersection square at a certain time = 806 465. It is characteristic 
of the stationary state; the particle configuration shows at certain points in 
space small fluctuating densifications which normally appear and disappear. 
However, the one in the circle starts acting as the nucleus of a jammed 
domain. Figs.[8)D-l8f show this domain at a succession of later times. While 
growing from the south and the west, it evaporates (but less fast) on the 
north and east side, with as a net result a displacement towards the south- 
west corner of the interaction square. Once it arrives there, it sticks to the 
entrance sites and blocks a set of horizontal and vertical inner lanes. We recall 
that the simulation includes the memory variables representing the lengths 
of the waiting line in each lane, even though these lines are not shown in the 
figures. In the course of time, the exact set of inner lanes that are blocked 
is subject to fluctuation but the jammed domain remains stable: we have 
continued the simulation until time t = 11 x 10^ without seeing it disappear. 
The final jammed state resembles closely a free flow state with an effective 
street width having a reduced value M' < M. 

The conclusion is that the continuous jamming transition observed in 
earlier work for M < 24 turns first order when M becomes larger. We 
have not been able to define a precise tricritical point (oc, Mc) at which the 
changeover takes place, but closer study, not reported here, reveals that the 
first order nature begins to set in as soon as M > 25. As a consequence, the 
straight line of data points in Fig. [7] cannot be continued beyond the range 
for which it is shown. 

5 Pattern formation and chevron effect 

Apart from the nucleation instability that it illustrates, Fig.[S] still shows an 
altogether different phenomenon that is of interest. It is the fact that in 
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Figure 8: Snapshots of the nucleation instabihty in a free flow state (see text). 
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the free flow regime the particles of the two types organize into alternating 
diagonal stripes [131 EH IIS]- Such stripes have indeed been observed in 
experiments and are also reproduced by realistic 'agent-based' models [HI 
\T7\ [TB] . In order to study them we replaced the particle dynamics with mean 
field equations in which on every lattice site r two continuous variables, 
p^{r,t) and p^{r,t), represent the densities of the eastward and northward 
traveling particles, respectively. These densities are postulated to satisfy 



where e^^y are basis vectors. As an auxiliary problem we solved these equa- 
tions on an intersection square with periodic boundary conditions, and found 
that the solution manifests the same stripe formation instability [13] as ob- 
served in the simulation of the true problem with open boundaries. We 
pointed out, however, that in the case of open boundaries the stripes in fact 
are not exactly at 45° with respect to the main axes, but that they form 
chevrons with a very weak opening angle, of the order of the degree. The 
effect appears again both in the particle simulation and in the numerical so- 
lution of the mean field equations (15.11) and may be observed in Fig. [91 The 
slope of the stripes has two roughly constant, but distinct, values the two 
triangular regions delimited by the dashed white lines. The chevron angle 
may be determined accurately by sufficient statistical averaging; it appears 
to be linear in a. 

In Ref. [13] we explained the origin of this 'chevron effect'. There is 
much to say about this phenomenon and a detailed account is in preparation 
[I4j : it includes cases with unequal injection rates in the two perpendicular 
directions, and with the intersection square subjected to cylindrical boundary 
conditions. In Ref. [TJ] we show how a particle may be localized by the wake 
of another particle of the same type, a mechanism which explains how global 
patterns are produced at the microscopic scale. 



Simple stylized models like the present one are certainly not meant to com- 
pete with more realistic ones [T71 [H]. The purpose of the model studied here 
is complementary. Since it may be analyzed more fully, it sets a standard 
scenario with respect to which others models may be discussed. Models of 
this type may also draw our attention to phenomena in traffic problems (of 
which the chevron effect is an example) that remain easily hidden in more 
elaborate many-parameter models. 

Real pedestrians in crossing flows have only two strategies available to 
avoid collisions, and experimental observation shows that they use both. 




[l-pf(r)]pf(r-e.)+pf(r + e.)pf(r) 
[l-pf(r)]pf(r-e,)+pf(r + e>f(r) 



(5.1) 



6 Outlook 
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Figure 9: Snapshot of an intersection square of linear size M = 640. The particles 
show stripe formation and a chevron effect (see text). Figure taken from |13] . 

The first one is adapting their speed; this is the strategy implemented in the 
present modeL The second one is deviating their trajectory; this would cor- 
respond to introducing the possibility of sideways motion. From observation 
we know that both strategies are used, but that large deviations from straight 
trajectories are relatively rare. We therefore consider the model presented 
here as a relevant starting point. Incorporating lateral motion is left for later 
work. 

Certain properties of the present model may not survive the introduction 
of sideways steps. In particular, the distinction between M successive jam- 
ming transitions in individual lanes is likely to get blurred; and if so, the 
question of the nature of the jamming transition will have to be asked anew. 
Other properties of this model, however, may well turn out to be robust. 
One example is the predominance of the flow through the outer lanes over 
those through the inner ones. Another one is the chevron effect; we believe 
that in future observations and experiments it will be worth looking for this 
effect. 
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